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ABSTRACT 
The  knowledge  of  a  complete  set  of  vector  modes  in  cylindrical 
waveguides  containing  inhomogeneous  media  not  necessarily  bounded  by  perfect 
conductors  permits  the  transformation  of  Maxwell' s  field  equations  into 
ordinary  transmission-line  equations.  The  determination  of  the  electromagnetic 
fields  produced  by  arbitrary  currents  in  such  regions  is  thereby  reduced  to  a 
conventional  network  problem. 
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1,     Introduction 

There  are  a  number  of  tX)imon  aspects  in  the  electromagnetic  problems 
that  arise  in  connection  with: 

a)  wave  propagation  along  dielectric  coated  surfaces  in  ■  open"  or 
■closed"  waveguides J 

b)  tropo spheric  or  duct  propagation  in  planar  (or  cylindrical) 
atmospheres^ 

c)  in  general,  wave  propagation  along  cylindrical  waveguides  that 
contain  inhomogeneous  media  and  are  not  necessarily  bounded  by 
perfect  conductors   '^  -• 

We  shall  be  concerned  with  the  determination  of  a  complete  set  of  vector  modes 
in  teniis  of  which  the  electromagnetic  fields  produced  by  prescribed  (or  induced) 
sources  in  such  inhomogeneous  regions  can  be  readily  evaluated.     The  desired 
modes  are  defined  by  a  vector  eigenvalue  problem  and  possess  orthogonality  pro- 
perties that  permit  the  reduction  of  a  general  source-excited  vector  field  pro- 
blem to  a  conventional  scalar  transmission-line  problem. 

In  the  steady  stete  of  time  dependence  e""^  ,  the  electric  field  E  and  the 
magnetic  field  H  produced  by  sources  of  electric  current  density  ^  and  magnetic 
current  density  M  obey  the  field  equations 

VxE     -     iiqiH     =»-M, 
(1.1a)  ^  - 

VxH    +    ikeE    -    J. 
where  k  «  co  ^/|J,  e     ,       In  the  cylindrical  regions  of  interest.  6  and  ti,  the  dielec- 

r       O     O 

trie  constant  and  permeability  relative  to  the  vacuum  values  e  and  a,   are  assumed 
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to  be  indepervlent  of  the  axial  coordinate  z  but  are  in  general  variable  over  the 
plane  transverse  to  the  z-axLs.  On  the  guide  walls  (if  present)  the  fields  are 
to  be  subject  to  the  boundary  conditions 


(1.1b)        E  -  oy»   H  X  V  , 

where  ^  is  an  impedance  dyadic  relating  the  electric  and  magnetic  field  components 
tangential  to  the  guide  walls.  As  noted  in  Fig.  1,  v  is  the  outer  unit  noimal  per- 
pendicular to  the  z  -axis  and  the  guide  walls  s^.  All  field  quantities  are  r.m.s. 
and  the  normalization  is  such  that  the  intrinsic  impedance  of  vacuum  is  unity  . 

The  explicit  solution  of  Eqs.  (l)  in  cylindrical  regions  may  be  effected 
by  a  representation  of  the  fields  in  terras  of  steady- state  modes.  The  latter  com- 
prise the  possible  guided  waves  that  can  be  propagated  along  the  axis  of  such  re- 
gions with  no  sources  present.  The  field  components  of  these  modes  are  expressible 
in  terms  of  vector  eigenfunctions  e  (p)  and  h  (p)  which  indicate  the 
(x,y)  5  p-dependence  of  the  source-free  electric  and  magnetic  fields  transverse 
to  the  transmission  direction  z.  Since  the  eigenfunctions  form  a  con5)lete  set, 
the  total  transverse  electric  and  magnetic  fields,  which  are  independent  field 
variables  in  Eqs.  (1),  may  be  represented  at  any  point  (x,y, z)  ■  r  as 

E.  (r)  -  TJS^)   e  (p) 
•^t  ^     *«-a  a    -~a  '-' 

(1.2) 

the  summation  being  extended  over  the  complete  set  of  modes.  With  a  proper  deter- 
mination of  the  mode  functions  e  and  h  the  original  vector  field  equations  (l) 

•~a  '-a  '^  ^ 

may  be  transformed  by  (2)   into  ordinary  scalar  differential  equations  for  the  mode 
amplitudes  V   (z)   and  I. (z): 

For  the  conventional  nonnalization  replace  H  and  ^  in  (l)  by 

-°     H     and    /5       J. 
^o    -"  /    «o       ~ 
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(1.3) 


dV, 

dT 

dl, 
dT 


ix  Z  I 
a  a  a 


ix  Y  V 
a  a  a 


-  V 


\' 


where  the  constant  parameters  x  ,  Z  •»  l/Y  and  the  z-dependent  source  functions 

V  ,  i  are  known.  Equations  (3)  pose  a  conventional  transmission- line  problem 
representative  of  Eqs.  (la)  with  the  transverse  spatial  dependence  of  the  typical 
mode  a'  suppressed.  This  problem  is  readily  solved  for  the  (voltage)  amplitude 

V  and  the  (current)  amplitude  I  ,  and  from  these  and  Eqs,  (2)  the  desired  repre- 

sentation  of  the  excited  electromagnetic  field  is  obtained.  It  should  be  noted, 

however,  that  the  ability  to  obtain  an  explicit  field  solution  is  predicated  upon 

the  explicit  evaluation  of  the  mode  functions  e  and  h  which  can  be  realized 

—a     --a 

analytically  only  in  cylindrical  regions  with  sufficiently  simple  e,  tx -variability. 

The  cylindrical  region  illustrated  in  cross-section  in  Fig.  1  represents 
the  general  uniform  waveguide;  the  transmission  direction  of  interest  is  along  z. 


Fig.  1.  Uniform  waveguide  with  variable  medium  e,  and  arbitrary  walls  (if  any) 


-  u  - 


A  special  case  to  be  treated  in  detail  is  the  planar  stratified  region 
depicted  in  Fig.  2a.  Although  pictured  as  a  region  within  which  the  dielectric 
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Fig.  2a.  Open  waveguide  with  planar 
stratification. 


Fig.  2b,  Closed  waveguide  with 
planar  stratification. 


constant  e  and  the  penneability  [j.  are  piecevn.se  constant,  both  e  and  ti  are  per- 
mitted to  have  a  relatively  arbitrary  variation  along  the  y-direction  of  stratif- 
ication.  A  region  of  this  type  is  frequently  termed  a  composite  open  waveguide. 
On  introduction  of  a  cartesian  coordinate  system, such  a  region  may  be  regarded  as 
a  uniform  waveguide  with  transmission  direction  along  either  the  x,  the  y  ,  or 
the  z  axis.  On  use  of  a  polar  coordinate  system  the  region  may  be  regarded  as 
a  unifonn  waveguide  with  transmission  along  the  y-axis,  as  a  non-unifoim  radial 
waveguide  with  transmission  in  the  p-direction,  or  as  an  angular  waveguide  with 
transmission  in  the  9^-direction.  If  the  region  is  enclosed  by  planar  walls,  as 
for  example  in  Fig.  2b,  it  becomes  a  closed  composite  waveguide  of  rectingular 
cross  section.  The  boundary  walls  need  not  be  asstuned  to  be  perfect  conductors 
and  will  be  regarded  in  general  as  '  impedance  walls'  in  a  mariner  defined  more 
precisely  below.  It  is  of  interest  to  investigate  the  circumstances  under  which 
the  open  region  of  Fig.  2a  is  a  limit  case  of  the  closed  region  of  Fig.  2b  as  the 
walls  recede  to  infinity  since  the  determination  of  a  complete  set  of  modes  for  a 
closed  region  is  a  relatively  elementary  problem. 
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The  proper  determination  of  a  complete  set  of  vector  modes  in  the  above 
regions  poses  a  vector  eigenvalue  problem  in  the  cross-section  transverse  to  the 
transmission  direction.  The  spectrum  of  modes  obtained  is  usually  discrete  in  the 
case  of  closed  regions,  as  in  Fig,  2b,  and  generally  both  discrete  and  continuous 
in  the  case  of  open  regions,  as  in  Fig»  2a.  The  discrete  modes  characterize  the 
natural  (free)  cross-sectional  resonances  of  a  region  in  that  they  are  the  possible 
source-free  field  solutions  each  satisfying  the  prescribed  boundary  conditions  in 
the  cross-section;  they  constitute  the  proper  eigenf unctions  of  the  spectrum.  On 
the  other  hand  the  continuous  modes,  which  characterize  the  radiative  part  of  the 
spectrum,  individually  are  possible  field  solutions,  but  only  as  a  complex  (a  wave- 
packet)  do  they  satisfy  the  required  cross-sectional  boundary  conditions;  they  con- 
stitute the  improper  eigenfunctions  of  the  spectrum.  Of  the  discrete  modes  encount- 
ered in  both  open  and  closed  regions  those  characterized  by  an  exponentially  decreas- 
ing behavior  outside  the  dielectric  regions  are  designated  as  surface  waves.  They 
are  fxorther  characterized  by  a  relative  insensitivity  to  the  presence  or  absence  of 
guide  walls  in  the  dielectric  vicinity,  i.e.  to  whether  the  guide  cross— section  is 
closed  or  open. 

The  discrete  and  continuous  vector  modes  possess  orthogonality  properties 
over  the  cross-section  that  are  determined  by  the  nature  of  the  dielectric  medium 
cind  the  guide  walls,  if  any.  In  the  hermit ean  case  with  non-dissipative  dielectric 
and  'reactive'  walls,  the  orthogonality  properties  imply  no  power  coupling  between 
modes,  i.e.  the  total  power  in  the  field  is  the  sum  of  the  powers  carried  by  the 
steady- state  modes.  In  the  non-hermitean  case  with  dissipative  dielectric  and 

'impedance'  walls,  there  is  power  coupling  between  modes  despite  the  orthogonality 
properties. 

In  this  study  our  concern  will  be  both  with  the  determination  of  the 
complete  set  of  modes  in  a  fairly  general  uniform  guide  with  variable  e,  [i,    and 
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with  the  concomittant  transformation  of  the  general  field  problem  therein  into  a 
transmission- line  problem.  In  this  report  the  eigenvalue  problem  for  the  modes 
will  be  stated  in  general.  Explicit  evaluations  of  the  mode  f mictions  will  be 
performed  for  planar  stratified  regions  wherein  the  vector  eigenvalue  problem  is 
separable  in  a  subsequent  report  (Part  II).  The  methods  employed  are  a  generaliza- 
tion and  extension  of  those  reported  previously'-  -' . 

2.  Eigenvalue  Problem  for  Modes 

It  is  of  interest  to  ascertain  whether  or  not  the  possible  field 
solutions  in  a  source-free  but  otherwise  arbitrary  urdfoim  guide  suffice,  when 
suitably  superposed,  to  represent  the  field  vrhen  arbitrary  sources  are  present* 
In  the  absence  of  sources  the  possible  field  solutions,  each  satisfying  prescribed 
boundary  conditions  on  the  guided  walls^  are  designated  as  proper  modes.  From  gen- 
eral symmetry  properties  of  a  xmifoim  guide  one  can  infer  about  such  modes  a  number 
of  characteristics  that  are  independent  of  the  specific  nature  of  the  guide  cross- 
section  and  hence  are  applicable  to  a  large  class  of  unifoim  waveguides. 

The  distinguishing  symmetry  property  of  a  source-free  cylindrical  guide 
is  its  invariance  to  translation  along  the  guide  axis.  The  independence  of  the 
properties  of  the  medium  and  walls  of  a  unifonn  guide  on  the  axial  coordinate  z 
implies  a  field  dependence  characteristic  of  the  translation  operator  -r  ^  ♦ 
The  eigenfunctions  of  the  translation  operator  have  the  form  e   .  Accordingly, 
the  fields  of  a  possible  node  in  a  tiniform  guide  will  likewise  have  a  separable 
e   -dependence, with  x,  the  mode  wave  number,  to  be  determined.  The  steady-state 
electric  and  magnetic  fields  of  a  tjrpical  mode  are  thus  given  to  vri.thin  an  ampli- 
tude factor  by 

E  (r)  -  ^(p)  e^^\ 
(2.1) 

H  (r)  -?/(p)  e^^ 
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which  represents  a  guided  wave  of  invariable  cross-sectional  form  traveling  in 
the  direction  of  increasing  z,   Fvirther  information  is  derived  if  the  guide  is 
also  invariant  to  a  reflection  in  ary  cross-sectional  plane.  This  reflection 
symmetry  implies  that  under  a  transformation  z  -^  -z,  the  transverse  (t)  and  lon- 
gitudinal (z)  components  of  the  mode  fields  transform  to  within  a  phase  factor  as 

(2.2) 

Ejp,z)  ->  -Ej_p,  -s),        H^(P>z)  ->  H^(p,  z)  , 

Thus  if  Eqs.  (l)  provide  a  mode  solution  of  the  source-free  field  equations,  a 
'reflected  mode'  with  wave  number  -  x  and  with  field  components  related  by  (2) 
to  those  in  Eqs.  (1)  likewise  provides  a  field  solution. 

The  vector  mode  functions  ^(p)  and  ^(p)  together  with  x  are  to  be 
so  determined  that  (l)  provides  a  source-free  solution  of  the  field  equations 
(1.1),  i.e., 

(2.2a)  J      ^^ 

V  X  ?/-  -  -ike  C  , 

where  e  =  e(p)  and  |j.  =  ti(p)  and  V  «  V.  +  ixz  }   z  denotes  the  unit  axial  vector 
and  V,  the  transverse  gradient  operator.  In  accordance  with  (l.lb)  the  discrete 
mode  fields  are  to  be  subject  for  all  z   to  the  boundary  condition 


(2.2b)         /-|.^Xv 


on  the  guide  walls  s.     In  uniform  guides  with  reflection  symmetry  the  impedance 
dyadic  /  is  diagonal  in  a  s,  z  basis. 

An  alternative  form  of  Eqs.  (2)  is  obtained  by  decomposition  into  compon- 
ents transverse  and  longitudinal  to  z  (cf.  ref.  [2]).  The  transverse  component  ob- 
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tained  by  vector  product  raviltiplication  of  Eqs.   (2a)  by  z     is 

"^ 

t   z     -^  t     "^-"t   '-o 
(2.3a) 

and  the  longitudinal  component  obtained  by  scalar  product  multiplication  of 
Eqs.  (2a)  by  z  is 

V .  .  z    %  <^.   "   -iku,  "/V" 
t    "O   -t      ^2 

(2.3b)  ^^  ^ 

V.  .  Zr.^  z^  =  -ike  6,  . 
t   ^  t  ^o         z 

On  elimination  of  the  longitudinal  components  in  (3a)  by  use  of  (3b)  one 
Obtains  as  the  defining  equation  for  the  transverse  components 

X(f^  =  k(txl,.i  Vj  V^)  'U^-^^ 

xj^  -k(el,  .i  V,-  V^)  .  V^t' 
k 

where' 1.    is   the  transverse  unit  dyadic.     In  view  of  Eqs.    (2b)   and   (3b),   the  dis- 
Crete  mode  fields  of  Eqs.    (Ua)  are  to  be  subject  to  the  boundary  conditions 


V.    •     z    /■  (^  ^  *•  iku  Z"~(s)   z    x&,    •  \> 
t        "o     i::  t  ^     ss        '^o     '- 1      '^ 

(2.ab; 

V^   •  >A    >^    z     =  ike  Z     (s)-V-.   ;<  z     .  v 
t'—t        ^o  zz— t'-o       — 

on  the  cniide  walls  s,   Z       and  Z       being  the   x-independent  elements  of  the 
°  '~       ss  zz 

diagonal  dyadic  "^    •     Equations   (2)  or  their  transverse  equivalent,  Eqs,    (h), 
constitute  a  vector  eigenvalue  problem  for  the  determination  both  of  the  possible 
mode  fields     <L^t    /T.    arxl  of  the  possible  wave  numbers  x  in  a  uniform  region  with 
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variable  e,  y^  medium.  The  conventional  eigenvalue  problem  for  the  modes  in  a 
homogeneous  cylindrical  guide  with  perfectly  conducting  walls  is  a  special  Ccise 
of  Eqs.  ik)   in  which  e  "  1  "  ix  and  ^  ■  0.  In  the  present  case  the  possible 
mode  vectors  C^.  and  /y.    need  not  be  perpendicular  as  in  the  special  case. 

The  reciprocity  properties  of  the  electromagnetic  field  permit  a  ready 
derivation  of  the  orthogonality  properties  of  the  mode  fields  (l)  •  Two  cases 
are  to  be  distinguished:  the  hermitean  and  the  non-hermitean.  In  the  herraitean 
case  the  dielectric  medium  and  guide  walls  are  non-dissipative,  whence  e,  m., 
and  i^  are  real.  In  this  case  it  can  be  immediately  inferred  from  the  constancy 
of  the  real  power  flow,  Ret^^      *  z  >  along  the  guide  axis  that  the  possible 
squared  wave  numbers  x  of  the  mode  fields  (1)  must  be  real.  If  two  such  modes 
are  distinguished  by  subscripts  a  and  b,  the  complex  reciprocity  theorem  for 
source-free  fields 

(2.5)  V  .  (E^xH^  +  E^)(H^)  -  0 

implies  that 

V  •  ^C^^'^l^g%7^)   +  i(x-x*)  z   .  iS^.^l*    ^!./^J  -0, 
X   r-a-^b   '^b'-a      ab  ^^   v-  a  ~  b    '^  b  —  a 


Cn  integration  over  the  guide  cross-section  S  and  use  of  the  two-dimensional  diver- 
gence theorem']  one  obtains,  in  view  of  Eqs.  (2b)  and  the  realness  of  the  diagonal 
dyadic  i  ^  » 

i(x-x*)  //((^  -^xz^  +  -^.'^n"  <^l^   dS 
a  o   y/  y       iia-^b  ~o    ^  a.    -*-  o  ~  b 

S 

■  -/K  X  V  •  (^,  -  I  •  ?/,  X  V)  ,  ?/,  V  V  -(^.-l.  V,  "     V)*| 

The  following  discussion  of  orthogonality  properties  is  similar  to  that  given  by 
Adlerl-  ^  and  is  included  because  of  its  relevance  in  the  anisotropic  cases  to  be 
discussed  in  a  separate  renort  (Part  III). 

/  /  V. 'A.  dS  ■  /y   •  A.ds  when  A  is  suitably  continuous  and  the  line  integral 
with  respect  to  s  is  taken  over  the  periphery  s  bounding  S. 


ds  "  0 
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whence  by  Eqs,    (2b)  one  has  the  bi-orthogonality  property 

(2.6a)  ffi^^K  X   z     -^  V   .   z  ^  <f^  J)  cJS  =  0       if     X     /    X*. 

S 

In  addition,  if  instead  of  the  mode  a  the  'reflected  mode'  ^  is  used  one 
obtains  by  (2)  the  orthogonality  property 


s 

On  addition  of  Eqs»  (6a)  and  (6b),  one  has  the  simple  bi-orthogonality  property 


S 

In  the  non-hermit ean  case  the  dielectric  medivun  and  guide  walls  may  be 
dissipative,  in  which  event  6,  p.,  X-      and  consequently  x  are  in  general  complex. 
The  orthogonality  properties  of  any  two  source-free  modes,  distinguished  by  sub- 
scripts a  and  b,  may  now  be  inferred  from  the  reciprocity  theorem 

(2.e)  V  .  ( E  X  H,   -  E.  X  H  )  =  0 

'^a   "^     '^   •**a 

in  a  manner  sijnilar  to  the  developnent  in  Eqs.    (5)-(7).     For  the  node  fields  defined 
by  Eqs,    (l)  one  derives 


//    ^  ^=*  ^K  '^  2«  -   f^'  ^^  <^J)  dS  =  0         if  -X     y  X 
v/  v/         ^  a     '^  b    ~vO        '^&    -o     ^  b  a 

S 

If  the  region  possesses  reflection  syminetry,  then 


S 


-  11  - 


whence  in  a   cylindrical  region  with  reflection  symmetry 


2    /     2 


(2.9)  .//       (f'.  y.xz     dS     -     0  ifxW 


X. 


a  '      b 


It  is  of  interest  to  observe  that  since  the  orthogonality  properties  (?)  and  (9) 
involve  only  transverse  field  components,  these  properties  may  equally  well  be 
derived  directly  from  the  transverse  eigenvalue  problem  of  Eqs,  (U), 

The  primary  importance  of  the  transverse  field  will  be  emphasized  by 
introducing  the  notation 

(2,10a)         jg^(p)  »  e(p)  and  ft^io)   =  T  h(p) 

for  the  transverse  components  of  a  mode  field,  the  parameter  Y  «  l/z  being  intro- 
duced to  secure  a  convenient  normalization  of  ^  relative  to  h.  On  substitution 
of  (10a)  the  transverse  eigenvalue  problem  (Ua)  can  be  restated  in  the  form 

X  Z  e  »  k(ti  1.  +  -iy  V.  i  V.  )  •  h  K  z 
^        t    ,c.       t  e   t   ^      -^o 

(2.10b) 

xY  h  -  k(e  1^  +  ^    V.  i  7^)  •  z  >^  e 

'-  t    ,2   tut    —o       '— 

subject  to  the  boundary  conditions  that  follow  from  (Ub)»  As  observed  in  (2), 
there  exist  both  +  x  and  -  x  eigensolutions  of  Eqs,  (lOb),   If  the  parameter  +Y 
is  associated  with  the  +  x  solution  and  -Y  with  the  -  x  solution,  then  it  is 
evident  from  either  (2)  or  (lOb)  that  the  mode  functions^  ji  for  the  ±   x  solu- 
tions become  identical.   In  view  of  this  identity  as  well  as  of  the  orthogon- 
ality properties  (7)  and  (9),  a  mode  will  be  distinguished  hereafter  not  by  the 
wave  number  x  but  rather  by  x''.  Such  a  mode  is  generally  a  standing  wave  formed 
by  superposition  of  a  wave  traveling  in  the  positive  z-direction  and  its  reflec- 
tion traveling  in  the  negative  z-direction.  The  transverse  field  components  of  a 
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2 

typical  mode  x  in  a  source-free  region  can  thus  be  represented  as 

^^(r)  =  V(z)  e(£) 

(2.11) 

H.(r)  »  I(z)  h(p) 

where  by  Eqs.    (l)   and  (2)   the  z-dependent  mode  amplitudes  V  and  I  have  the  fonn 

,   .  ixz  ^   ,     -ixz 

V(z)  -  A^.  e         +  A_e 

I(z)  =  Y(A_^  e         -  A_e         )  e 

In  consequence  of  the  arbitrariness  of  A  +  there  is  an  arbitrariness  in  the  mode 
amplitudes  V,  I,  which  is  exhibited  in  a  significant  form  by  noting  that  V  and  I 
satisfy  the  first-order  differential  equations 


(2.12) 


dz 


~  -  ixY7, 
dz 


where  Z  =  l/T.  Since  Eqs.  (12)  have  the  foim  of  homogeneous  transmission-line 
equations,  one  employs  the  conventional  designations  of  V  as  the  'voltage'  ampli- 
tude and  I  fs  the  'current'  amplitude  of  the  mode  in  question. 

In  an  arbitrary  cylindrical  waveguide  excited  by  sources,  the  ability 
to  represent  the  electromagnetic  fields  as  a  superposition  of  source-free  modes 
of  the  form  (11)  is  predicated  upon  the  completeness  of  the  set  of  possible  mode 

functions  e  ,  h  .  The  completeness  will  be  exhibited  in  the  subsequent  procedure 
~*a'  *^ 

employed  to  obtain  explicit  solutions  of  the  vector  eigenvalue  problem  (cf.  Part  II)« 
The  development  to  this  point  has  indicated  that  if  in  a  general  source-free  region 
possessing  translation  and  reflection  symmetry  there  exist  possible  mode  functions 
Bj   Jj.  defined  by  Eqs.  (10b),  they  must  possess  orthogonality  properties  of  the 
form  (7)  or  (?)• 
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3»     Field  Representations  in  Regions  vrith  Sources 

The  vector  modes  characteristic  of  a  source-free  region  provide  a 
represent? tion  of  the  electromagnetic  field  in  the  pi^sence  of  sources.     In  this 
section  we  shall  describe  such  a  representation  for  a  unifoim  wa^veguide  with  arbit- 
rary cross-section,  with  arbitrary  isotropic  e,  \i  cross-sectional  variation,   and 
with  rather  general  boundary  conditions  at  the  guide  walls   (if  any).     The  electro- 
magnetic fields  are  determined  by  the  inhomogeneous  field  Eqs.    (1,1a)  subject  to 
the  boundary  conditions  (l,lb)  on  the  guide  walls.     As  in  Eqs. (2»2)  -   (2.U)  one 
decomposes  the  inhomogeneous  field  equations  into  transverse  components 

V^  E     -  4-     E .    =  -iku  H.  X  z     +  M  ><  z 
t     z       3z    -t  ^  <-t    -o      .-    *-o 

(3.1)  g 

and  longitudinal  components 

V.  •  z  X  E.  -  -iku  H  +  M 
t   '-^   ~^t     ^  z    z 

(3.2) 

t   --H   «-o        z    z  * 

whence  on  elimination  of  the  longitudinal  comporvants  from  (l)  one  obtains 
as  the  defining  equations  for  the  transverse  field  components: 

1-  E,  =  ik     f p.  1.   +  i    V.  i  V.  A    •  H.  X  z       +    M  -A  z     +  iV.   (  T-r-J 


(3.3) 


d 
S7 


H.    -  ik    f  e  1^  +  i     V.^vA    *   z  xK     +  z    '^    J  +  i^+    (  TT 
-'t  I        t       .  2       t  n     t  y       --0    — t        -o      '^  t   y  )s\ilJ  > 


subject  to  the  bovuidary  conditions   (l»lb)» 

Since  the  vector  modes  defined  in   Section  2    can  be  shown  to  comprise 
a  complete  set,   they  siiffice  to  describe  the  electromagnetic  field  in  a  guide 
containing  sources,     As  the  characteristic  field  patterns  e  (p),  h  (p)  of  each 


-  in  - 


mode  are  known  over  the  guide  section,  only  the  mode  amplitude  variation  along 

the  guide  has  to  be  determined.  Accordingly,  the  transverse  r.m.s.  components 

of  the  electric  and  magnetic  fields  will  be  represented  by  the  linear  superposition 

(3.U) 

a 

here  and  in  the  following  the  summation  is  to  be  extended  over  all  possible 
values  of  the  mode  subscript  a  in  the  discrete  and  continuous  spectinun. 
a)  Hermitean  case 

In  the  hermitean  case  wherein  the  guide  has  a  non-dissipative 
median  and  reactive  walls,  the  mode  functions  e  and  h  can  be  so  normalized 
tliat  their  bi-orthgonality  properties  can  be  xjritten  in  accord  with  Eq,  (2,7)  as 

0,    a     ^    h 


(3.5)        //  e  .  h*v  z   dS  «  5  ,   «= 

yy      ~a  -^    --o  ^^    1    a  -  b, 


the  surface  integrals  being  extended  over  the  entire  guide  cross-section  S,     Eqs, 
(5)  imply  that  the  total  complex  power  flow  at  any  point  z  is  given  by 


(3.6) 


P     =     /  /     E.    X   hJ  .     z^     dS     "    Y    V  (z)   I*(z) 
yj       ~t       ~t        -0  ^       a  a 


and  is  the  sim  of  the  powers  carried  by  the  individual  modes.  In  view  of  the 
orthogonality  properties  (5)  the  mode  amplitudes  and  the  transverse  field  are 
related  by  the  transform  equations 

s 

(3.7) 


v^^  "// v^^  '^o^4^^^ 


65     • 
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The  explicit  determination  of  the  mode  amplitudes  V  ,  I  is  accom- 

plished  by  transformation  of  the  transverse  field  equation  (3)  in  accordance  with 
the  operations  indicated  in  Eq.  (?)•  Thus^  integrating  over  the  guide  cross- 
section  the  scalar  product  of  the  first  of  Eos,  (3)  with  h  x  z  and  the  second 

<~a    '~o 


with  z   X  e  ,  one  finds  by  Eqs »  (7)   (cf.   ref , H ) 


'o    '^a. 


d7 

a 


dz 
(3.6e)         ^j 


ik  //  h*x  z     •   (n  1.  +  i    V.   i  7.)    »  H.  >  z     dS  +  v  (z) 
v/  J     ~-a    -^xD         ^    t     ,  2       t  e     t        '-'t    '"o  a 


where 


(3.8b) 


dT  -     ^//  -^o^^a  •   ^^  V  ^     \l  \)    •  ^^  5t  <5S  *  ijz) 
t 

S  S 


s 

The  integrals  on  the  right-hand  side  of  Eqs.    (8a)  may  be  simplified  on  use  of 
integration  by  parts,   which  gives  the  adjointness  relation 

//       (^^  -  ^^   •  H^><  z  )   •   i^*'^  z     -H^xz     .    (V^  i  7.«  h*xz   )     dS 
J  J      [_    t  e     t       ~t    --0         ~a     -^0       -^t     -o  t  6     t     -^a    ~o  J 

S 

-     /    -     (v  •   h*  X  z   )  V.  .  H.  X  z     -  (v  .  H.  ><  z   )  7.    •  h*  X  z       dS 
/      e    [_  ~      -a    -^o'     t    '-'t     "O         ~       ~t    -^o       t      ^a    '^oj 

s 

and  the  dual  relation  obtained  by  the  replacement  of  H. ,  h  ,  e  by  E. ,  e^,  n, 
respectively;  the  surface  integrals  over  the  guide  cross-section  S  are  expressed 
in  terms  of  line  integrals  over  the  peripheral  curve  s  bounding  the  cross-sectioni 

'  The  quantity  i  is  of  course  not  to  be  confused  with  the  symbol  i  for-/3. 
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In  view  of  the  defining  Eqs,  (2.10b)  for  the  mode  functions  e  ,  h  and  the  bound- 
ary  conditions  (1,1b),  the  line  integrals  in  Eqs.  (9)  vanish  and  the  transformed 
field  equations  (8a)  become  by  Eqs.  (?) 


dV 

--^  »  ix  Z  I  -  V 
dz       a  a  a    a 

(3.10) 

dl 

:—  =  ix  T  V  -  i  . 
dz       a  a  a    a  * 

The  introduction  of  the  mode  representation  (1|)  has  trsnsformed  the 

Original  vector  field  problem  posed  by  Eqs.  (l.la),  (1.1b)  into  the  •  a- dependent ' 

scalar  differential  equations  (10),  The  latter  pre  conventional  inhomogeneous 

transmission-line  equations  of  the  same  nature  as  those  encountered  in  the  mode 

treatment  of  homogeneous  perfectly  conducting  unifoim  guides.  The  determination 

j^  the  behavior  of  the  voltage  V„  and  current  I^  is  a  straightforward  network  pro- 

blem  of  finding  the  response  to  the  known  'source  voltage',  v^,  and  'source  current', 

i  .  With  the  evaluation  of  V  and  I^  from  Eqs.  (10)  the  electromagnetic  field  at 
a  a     a 

any  point  r_  can  be  synthesized  by  vse   of  Eqs,  (li). 

The  knowledge  of  the  dependence  of  the  x  ,  Z  ,  v  ,  and  i  on  the  struc- 

51     cl     cL  a 

ture  of  the  uniform  guide,  and  on  the  frequency,  the  excitation,  etc. ^ is  a  necess- 
ary prerequisite  to  the  solution  of  Eqs.  (10),  The  quantities  x  and  Z  are  deter- 
mined  from  the  solution  of  the  eigenvalue  problem.  The  dependence  of  v  and  i  on 
the  excitation  J  and  M  can  be  cast  in  a  form  more  general  than  that  of  Eqs,  (8b) 

since  the  latter  assumes  differentiability  of  J  /&   and  M  /n  over  the  guide  cross- 

z       z 

section.  The  desired  form  is  obtained  on  integrating  the  right-hand  integrals  in 
(8b)  by  parts  and  assuming  the  J„  vanishes  on  the  guide  periphery;  one  finds 
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^=(z)     =     //^^M    •  h*(p)   dS     +     Z     //  J(r)    •  e*  (p)  6S 
S  S 

(3.n) 


a 


(z)  =  ff  jjs)  •  «:(p)  ^  -  \//^^^^  •  4(5)  <^. 


where  by  definition 


V. -h     X  z 

Z     e         =     -  ■■'■  ■  ■.■; —  z 

a  '^za  IKS         -^o 

V,»z    X  e 
V     V,         =  t  -^o      ^a 

in         =     -  .X  z  . 

a  ^za  zky.         --o 

The  expressions  (ll)   for  v     and  i  ,  which  sre  applicable  even  for  discontinuous 
J  /e  and  M_/|J-,   are  fundamental  to  the  netv/ork  calculation  of  the  fields  produced 
by  arbitraiy  excitation  J  and  M   in  an  arbitrary  cylindrical  guide, 
b)     Non-He CTiite an  case. 

When  the  medium  and  guide  walls  are  dissipative,   the  vector  mode 
functicnse     and  h     will  be  so  normalized  that  their  bi-orthogonality  properties 

a,  S 

may  be  written  as 


(3.12)  //    g^  .  ^  X  ^^  dS  -  6^^ 


in  acccord  with  Eqs.  (2,9).   In  contrast  to  the  hermitean  case  the  complex  power 
flow  cannot  be  expressed  as  the  sum  of  the  complex  powers  carried  by  the  individual 
modes,  Ne-rertheless  the  treatment  of  this  case  parallels  that  presented  in  Eqs, 
(5)- (10)  for  the  hermitean  case  provided  only  that  the  complex  conjugate, desig- 
nated by  the  asterisk,  is  omitted  in  these  equations.  With  this  omission,  the 
transmission-line  equations  (10)  are  still  applicable  in  this  case  except  that 
v  and  i  are  determined  by 
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V  . 
a 


^(z)  "//^{r)   .  4(p)  dS  -  Z^  //  J(r)  .  e^^(p)  dS  , 
S  S 

(3.13) 


1 
a 

S  S 

which  differ  from  (11)  only  in  that  the  sign  of  the  last  term  is  negative,  since 

e  and  h  are  imaginary, 
z      z 

li.  Planar  Stratified  Regions 

The  deteimination  of  the  modes  in  an  arbitrary  uniform  waveguide  with 
axis  along  the  z-direction  and  with  variable  e  and  |i  in  the  cross-section  trans- 
verse thereto  leads  to  a  vector  eigenvalue  problem  which  is  in  general  difficult 
to  solve.  For  the  special  case  of  a  planar  stratified  cross-section  (cf.  Fig,  2), 
wherein  e  and  (x  are  scalar  functions  of  only  one  transverse  coordinate  y,  the 
vector  eigenvalue  problem  can  be  separated  into  two  scalar  eigenvalue  problems. 
This  separability  is  not  \inexpected,  for  such  a  region  can  be  regarded  not  only 
as  a  guide  with  transmission  axis  along  z  but  also,  for  the  purpose  of  mode  de- 
teimination, as  a  simple  teiminated  unifoim  guide  with  axis  along  y.  When  the 
functional  fom  of  e(y)  and  n(y)  is  sufficiently  simple,  explicit  expressions  for 
the  desired  mode  functions  can  be  obtained  by  evaluation  of  the  transverse  reso- 
nances in  this  terminated  waveguide.  The  latter  are  defined  by  a  general  Stuim- 
Liouville  problem  that  admits  both  a  discrete  and  continuous  spectrum  and  will  be 
evaluated  explicitly  in  another  report  (Part  II)  for  a  number  of  typical  cases. 

The  vector  eigenvalue  problem  defined  in  Eqs,  (2,2)  can  be  decomposed 
into  components  transverse  and  longitudinal  to  the  y-direction  of  stratification. 

Then  by  a  development  similar  to  that  of  Eqs,  (2,3)-(2.U)  one  obtains  as  vector 
eigenvalue  equations  for  the  desired  mode  fields  transverse  (T)  to  the  y-axis, 
that  is,  for   ^^(Xjy)  •  '^'  and  ^^(x,y)  e  '^, 


(U.la) 
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^^  /  1  X  ^ 

^    -    die  (eV    -r    Vt^  •   Jo^-^T 


with  the  longitudinal  components    £     and    j?y     given  by 
-  ^  ^y  -     \  •  Jo^  #T 

where  V„-V-y      »-    -     x      ^    +     !'<■  z     and  e  =  e(y),     n  "  l^(y)«     In  accord 

with  Eqs,    (2.2b)   the  discrete  mode  fields    ^  and    /^  are  to  be  subject  to  the 
boundary  conditions 

(k.2a)  /^»+zV-         and<!^»     +     Z         Pt 

^"•^    '  ^y  —      yyz  z  zzy 

on  the  x-boxindaries   (if  ary)  of  the  cross-section,   and  to 

(U.2b)  4     -"    2^     ^z       ^      ^z     -    1    ^zz-^x 

on  the  y-boundaries  (if  any).  The  impedance  coefficients  Z^  are  prescribed 

constants  but  may  be  different  on  the  x-  and  y-boundaries;  the  +  signs  refer 

respectively  to  conditions  at  the  boundary  with  the  larger  or  smaller  x-  or 

y-coordinate. 

For  the  special  esse  wherein  the  cross-section  is  \inbounded  in  the  x- 
direction,  there  exists  a  set  of  modes  independent  of  x  that  suffice  for  field 
problems  in  v*iich  d/dx  -  0.  Under  such  circumstances  the  vector  equations  (l) 
automatically  decompose  into  two  distinct  and  separate  types.  The  mode  field 
of  one  type  are  found  to  be 
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k> 


or  on  elimination  of   fh  , 
'  z  ' 


(U.3b)  (9      1      I-    +     k^e-i^  )        (f-0, 

where  by  (2b) 

on  the  y-boundaries.      It  is  to  bo  observed  from  Eqs.    (lb)   that  k|i, /V    «     <^,and  that 

y        X 

hence  Eqs«  (3)  characterize  'H-type'  modes  that  possess  an  H-  but  no  E-  field 
component  along  the  y-direction.  The  mode  fields  of  the  second  type  are  con^s- 
pondingly  determined  by 

(u.ba)  ^  -  iJc(^L  -  -|-  )  ;V^^ 


or,  on  elimination  of    O  ,  by 

Li 

(u.Ub)  (3    1  |.  +  k^^.  2^  )  ;^.o, 

3y     e     dy  ^         e        '    x         » 

where  by  (2b) 

{\u\xc)  £    -  +    Z       7^ 

Z         —         ZZ  '   X 
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on  the  y-bound  ries.     From  Eqs.    (lb)  one  notes  that  -ke  (S    ■  X/^  and  that  hence 

y  * 

Eqs.    (U)  characterize    'E-type'   modes  with  an  E-  but  no  H-  field  ccmponent  along 
the  y-direction,     As  in  the  previous  type,    this  designation  is  appropriate  to 
the  view  of  the   cross-section  as  a  terminated  waveguide  with  transmission  along 
y.     From  Eqs,    (2,7)   and  (2,9)   the  possible  mode  solutions    S  ,1^  oi  Eqs.    (3) 
and   (Ii)   are  known  in  the  hermitean  case   to  possess  orthogonality  properties  of 
the  form 

S 

where  the  integral  is  extended  over  the  original  xy  guide  cross-section  and  where 

(«) 

the  conjugate  symbol         is  to  be  omitted  in  the  non-hermitean  case. 

I'Jhen  there  is  field  variability  along  x,    the  vector  eigenvalue  problem 
(l)  does  not  automatically  separate  into  the  two  distinct  types  given  in  Eqs»   (3) 
and  (U).     However,   since  the  stratified  region  can  be  regarded  as  a  uniform  wave- 
guide with  homogeneous  cross-section  and  with  axis  along  y,   the  conventional  mode 

/^  ix  s 

formalism  can  be  employed  to  reduce   the  determination  of  the  mode  fields   ^e 

and    ^e  to  a  transmission-line  problem  for  their  essential  variation  along  y. 

With  respect  to  the  y-axis  one  can  express  the  transverse  field  components  of  a 
typical  mode  as 


^Jr)     -    V(y)     £(x,z) 
(U.6a)  ^  ^ 

A^^(r)     -     i(y)     h(x,z) 


whence  by  Eqs.   (lb)   the  longitudinal  components  follow  as 
(U,6b) 


-iJce    <^y     -     I(y)     ^T  •  ^  '^  ^o 


-iJqi   Vy     -     V(y)     V^  •  y^  '^    e   . 
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The  raode  frnictlons  e,  h  are  to  be  so  determined  as  to  satisfy  the  specified  con- 
ditions  (2a)  on  the  x-boimdaries,  and  the  mode  amplitudes  V,  I  are  to  satisfy 
the  conditions  (2b)  on  the  y-boundaries. 

For  the  special  case  of  perfectly  conducting  x-boundaries ,  Z   ■  0  ■  Z 
in  (2a)  and  the  mode  functions  e,  h  separate  into  the  conventional  TM  (E-)  and 
TE  (H^)  types  with  respect  to  the  y-axis.  The  TM  modes  are  defined  by  (cf,  ref,  [2]) 


(U.7a)  a  •(x,z)     -  -V„  gfix)  e^^^  =  h'x  y^ 

with 

(U.7b)  (d        +     k'    )  <jf{x)  =  0 

dx^  ^ 

subject  to  ^  -  0  on  the  x-boundaries.  The  TE  modes  are  defined  by 


(i+.8a)  h"  (x,2)  =  -V„  T(x)  e^""^  -  y  y  S" 


A.0 


with 


2        2 


(U.8b)  (^  +  k   )  Y(x)  =  0 

dx      ^ 

subject  to  dl/dx  =  0  on  the  x-boundaries.     On  substitution  of  Eqs,    (6)-(8)  into 
(la), the  amplitudes  V,   I  for  both  mode  types  ar«  found  to  obey  the  transmission - 
line  equations 


(U.9) 


^    -     ix  Y  V 
dy 


where  the  variable  propagation  wave  number  x  and  characteristic  impedance  Z  »  1/T 
are  given  for  the  1M  modes  by 
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(U.lOa) 


Z 


X 


X 


/ — 2 2 ? 

y  k  eti  -  k^  -  X  J 


and  for  the  TE  modes  by 


(U.lOb) 


<v"       1  k 

Y  X 


(note  k^-  k^  •  k^). 


"j/     kcM.-k     -X 


In  view  of  (6a)   the  mode  amplitudes  V,   I  are  to  satisfy  in  addition  to 
Eqs,   (9)   the  conditions   (2b)  on  the  y-boundaries»     For  the  case  wherein   a/  is 
constant  on  the  y-boumdaries  and  the  x-boundaries  are  perfect  conductors,  the 
boundary  conditions   (2b)   can  be  phrased  in  teniis  of  a  combination  of  IK   and  TE 
modes,  viz: 


(U.lla) 


/v|     ^t  /Ml       ^H 

7    e      +    V      e 


,/M  I  Aft  «     ^ 

(I    e      +1      e     ) 


on  the  y-boundaries,or,  equivalently. 


(U.llb) 


t  II 


t  ti 

Z21I      +     Z22     I 


where 


^■y7(i'-f-s''^'  ^2-//(^~''>-"'^-^2i.  ^.3=//(r-/-i")as. 


A  further  specialization  arises  in  the  case  of  homogeneous  isotropic  y-boundariesj 
in  this  event  the  coupling  of  the  TM  and  TE  modes  implied  by  (11)  is  removed  since 

\2  -"  '21  -  ^  -^  ^11  -  ^22- 

The  transmission-line  equations  (9)   subject  to    (ll)   constitute  a 

•transverse  resonance'    problem  in  the  xy-cross-section.     Mathematically,   they 
characterize  a  one -dimensional  eigenvalue  pitJblem  of  the  Stuwi-Liouville  type. 
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For  the  case  of  the  TM  modes,  elimination  of  v  from  (9)  leads  by  use  of  (lOa) 
to  the  scalar  eigenvalue  equation 


y^W 


(U.12a) 


(^U*>=V-^)   ^^'w  ■°- 


/^ii 


Elimination  of  I       from  (9)  and  use  of  (3D)  leads  to  the  scolsr  TE  eigenvalue 

equations 

2     2 
(U.12b)  ^    i    f    +     k^e-     JlL)    v"(y)  -0. 

For  the  case  of  imperfectly  conducting  and  non-isotropic  y-boundaries  the  eigen- 
value equations  (12a)  and  (l2b)  are  coupled  by  the  boundary  conditions  (ll).  As 
noted  above,  this  coupling  disappears  for  isotropic  walls  and  for  perfectly  con- 
ducting walls.  In  the  latter  event  the  boundary  conditions  on  Eqs.  (12)  are, 
respectively. 

At 

(u.i3a)  It  "  ° 

and 

(h.l3b)  v"  "  0 

on  the  y-boundaries.  One  observes  that  the  eigenvalue  problems  (l2a)  and  (12b) 

encountered  in  the  case  •sf—  ■  0  are  not  essentially  different  from  those  encoun- 

dx 

tered  in  the  cases   (3b)   and   (Ub),   where  d/dx  "  0. 

The  detennination  of  the  solutions  to  the  eigenv.^lue  problems  posed  by 
Eqs.    (7b),    (8b)   and  Eqs.    (12a),    (12b)  may  be  accomplished  by  the  resolvent  or 
characteristic  Green's  function  method   L  J ,     Since  the  eigenvalue  equations    (7b), 
(8b)  are  special  cases  of  Eqs.    (12),   only  the  latter  will  be  considered.       The 
characteristic  Green's  function  g(y,y  ;   X)   associated  with  Eqs.    (12a),  for  example, 
is  defined  by 
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u.u)  ^h  ^  h  *  ^^^''  ¥  g(yV>>^) --6(y-y'). 


Let  fLCy)  =  ^a^y)  denote  the  desired  eigenfunctions  normalized  with  weight  func- 

2         2 
tion  l/e(y)  to  unity  over  the  y-interval,  and  let  X    "  ^v  *  ^a      denote  the  desired 

eigenvalues.     Then  (cf»  ref,  [l] )  one  has 


(4.15) 


;(y').,5(y-y')  -  -^    y^    g(y,y'j  x)  d\    -   ^I   ^a^?)  ^(y), 


where  the  contour  integral  is  to  be  extended  about  all  singularities  X     of  g  in 
the  complex  X-plane  and  where  the  s\im  over  a  includes  contributions  from  the  com- 
plete discrete  and  continuous  spectrum.     As  before  the  complex  conjugate  symbol 
in  Eq,   (l5)  is  omitted  in  the  non-hermitean  case.     Since  the  characteristic  Green's 
functions  defined  by  the  inhomogeneous  differential  equation  (lii)  plus  appropriate 
boundary  conditions  can  be  explicitly  deteimined  if  the  functional  dependence  of 
e(y)   and  |j.(y)  is  sufficiently  simple,  Eqs,    (l5)   constitute  a  pragmatic  and  formal 
solution  of  the   scalar  eigenvalue  problem  as  will  be  illustrated  in  another  report 
(Part  II,) 
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